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AN EMPIRICAL FORMULATION OF THERMIONIC CONVERTER PERFORMANCE 
AS A FUNCTION OF ELECTRODE EMISSION PROPERTIES 
by John R. Smith and Arthur L. Smith 
Lewis Research Center 

SUMMARY 

A simple empirical equation giving the emission current density in terms of the 
cesium atom arrival rate and the electrode surface temperature is presented in this re- 
port. The Swanson, Strayer correlation for the maximum change in work function due to 
cesium adsorption is compared with field-emission data for some refractory metals. 
Oxygen effects are discussed. 


INTRODUCTION 

Thermionic converter analysis would be facilitated by a simple equation giving the 
electron emission current density J explicitly in terms of the electrode temperature T 
and the cesium particle arrival rate p. 

The results should, of course, agree with existing experimental data. Also the ex- 
pression must predict J where no experimental data exist (e. g. , for p larger than in 
existing data and for untried electrode materials). 

Even with the progress that has been made toward developing a general understand- 
ing of the effect of alkali metal adsorption on the electron work function, there is as yet 
no basic calculation of J as a function of T and p. A number of semiempirical cor- 
relations (refs. lto8), however, relate J to T and cesium gas properties. Usually 
the most tested and discussed of these are the classical atomic physics approach of 
Rasor and Warner (ref. 3) and the molecular chemistry formulation of Levine and 
Gyftopoulos (ref. 5). 

Although both works have in some instances agreed with experiments, some basic 
assumptions of the Rasor, Warner and the Levine, Gyftopoulos correlations have been 
questioned (refs. 9 to 12), both disagree with some experiments (refs. 13 to 18). While 
these works at least can be used as complex data fits, straightforward empiricism can 
yield a simpler formulation. 



Presented in this report is a simple, accurate converter physics formulation. The 
result is an explicit equation for J in terms of the bare work function cp q, p, and T. 
Since the derivation is purely empirical, agreement with experiment is improved. This 
formulation is used to analyze the performance of an oxygenated cesiated diode. 


THE CATHODE 

A good deal of information has been obtained empirically about converter cathode 
surfaces. For example, Breitwieser's experiments (ref. 19) yield nearly linear plots 
of In J against In Pq s , cesium pressure, for constant T’s. In addition, the plots of 
effective work function cp against T/T r are nearly independent of T f for experiments 
with p = P Cs / \/27rmkT, where P Cs is the vapor pressure of cesium at T r - And cp, 
T/T r relations are linear for 2. 3 < (p < 3. 3, the useful range for power producing 
thermionic converters. Further, experimenters (ref. 20) provided methods to reduce 
data scatter below that of 0. 1 eV for cp, T/T r plots . 

Houston (ref. 17) using the ’’plasma anode” technique, has shown that each of the 
metals tungsten (W), palladium (Pd), ruthenium (Ru), osmium (Os), rhodium (Rh), 
iridium (Ir), and platinum (Pt) gives its own curve when cp is plotted against T/T r for 
cesium (Cs) reservoir temperatures of 414, 434, and 454 K. These data are nearly 
linear for 2. 3 < <p < 3. 3 eV and can be correlated by 

<?= — + b (1) 

T r 

[-(aT/T +b)kfl 

J = 120 T exp 1 J (2) 

where a and b are to be determined empirically. Equations (1) and (2) apply where 
p = P Cq /v/2St;. Generalizations will be discussed in the section THE CATHODE. 

If experimental data were available for all possible cathode surfaces, it would be 
merely necessary to fit for a and b. But, of course, such data are not available; so it 
is desirable to write empirically a and also b in terms of cathode properties. The 
bare work function cp^ is a property which distinguishes cathode surfaces; so it is de- 
sirable to find a (cpq) and b(<pQ). 

Since the work function of a surface is very sensitive to surface condition such as 
contamination, history of heat treatment, sputtering, and bulk impurity content, the fol- 
lowing criteria should be satisfied by experiments designed to determine a (cp^) and 
b(o? Q ) empirically. First, it is desirable to have all samples measured in the same tube. 
Second, the bare work function should be obtained in the same tube that the cesiated data 
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are obtained. Houston (ref. 17) and Wilson (refs. 21 and 22) obtained experimental data 

in this manner. Their results are not in agreement, however. Houston offers at least 

plausible arguments for some disagreement in Wilson’s results; so Houston’s data are 

used in this report. The values of a and b determined empirically can, of course, be 

improved as progress is made experimentally. 

A least-mean-square straight-line fit was used to determine a and b for the 

metals reported by Houston (except for which there were insufficient data). The results 

are given in table I. The rms deviations of aT/T + b from the data A____ and the 

r i ills 

magnitudes of the maximum deviations A mo „ are seen to be quite small. 

The relations developed are limited to poly crystalline surfaces of the metals studied 
in Houston's experiment. These developments may be made more general by taking note 
of Houston's plots of T/T r for a given <p plotted against <pQ. Inasmuch as these plots 
are nearly linear, an 1ms fit will produce a and /3, where 

— 0! (j 0r\ + (3) 

T r 


Extending this procedure to obtain two values of cp, say and cp 2 , would result in 


<Pl = + £j) + b 

> 

<p2 = a («2 Vo ^ ^ 


(4) 


And the solution for a and b as a function of cp^ would follow immediately. 

Good fits can be obtained by using equation (3) for essentially all values of cp within 
the range ~2. 20 to ~3. 4 eV (range of interest in all practical power producing thermionic 
converters). Plots of the data including the Rasor, Warner theory and the 1ms linear fit 
for .3. 3 and 2. 3 eV are shown in figures 1 and 2, respectively. The limitation of the 
method was determined by extrapolating the 1ms fit of the Pt group data down to = 

3. 5 eV and comparing these data with Houston's (ref. 17) results for thorium (Th) and 
hafnium (Hf), which were taken in a separate experiment using the "plasma-anode” 
technique. Although cathodes having bare work functions of less than 4. 3 eV may not be 
of interest for thermionic converters, it is interesting to note that at cp = 3. 3 eV the 
extrapolation agrees well with data on both Hf and Th, provided that extensions of the 
linear portions of the cp, T/T r curves are used. At 2. 3 eV, the agreement is good for 
Hf but not Th. The Rasor, Warner theory shows good agreement with the data at 3. 3 eV; 
however, at 2. 3 eV the Rasor, Warner theory does not correlate the data as well as it did 
at the higher work function. 
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If we let cp = 3. 3 eV and <p^ = 2.3 eV, equations (4) can be solved to yield 


a = 


0. 176 cp Q - 0. 158 
b = 3. 3 - a(0. 702 cp Q + 0. 193) 


(5) 




From the preceding, it appears that these results would be useful for 2. 3 < cp ^ 3. 3 eV 
and perhaps for 4. 3 < <pQ < 6 . 0 eV for Cs on metal substrates. 

As a test of the self-consistency of this procedure, T/T r = &%(pQ + was fitted at 
= 2. 8 eV. Then ac^a^cpQ + ^ 3 ) + b was compared with 2. 8 eV and found to agree to 
within ±0. 05 eV for 4. 1 < cp^ < 6 . 5 eV. Table n gives the deviations A, s using equa- 
tions (5). They are larger than in table I, but the A 's are all less than 0. 1 eV. 

Equations (1) and (2) as written are not useful for converter system design. First, 
for a given T, T , and cathode surface, cp (and hence J) will generally depend on the 
converter configuration, that is, the interelectrode spacing and wall temperatures. 
Second, it would be desirable to use the empiricism in cases where the adsorbate is not 
coming solely from a vapor, for example, where adsorbate is supplied from the sub- 
strate for impregnated cathodes. 

Close examination reveals that equations (1) and (2) could be ambiguous ( cp depends 
on more than T and T ) if it is assumed that cp and J are unique functions of T and 
P for a given cathode surface, where p is the cesium arrival rate. If, for example, 
an electrode is immersed in a gas of temperature T , pressure P^, and particle mass 
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m 


g’ 


P = 


/ / 0 » 

= P_/i/27rm kT . If Knudsen flow conditions apply, p = p n = evaporation rate at 

o r & & v — 

the cesium reservoir. But if continuous flow conditions apply, p = -/T r /T Pq. So for 
a given T and T , different p's and hence different J's are possible. 

In order to develop equations (1) and (2) in a generalized form (and consequently to 
make the usefulness of a's and b's more general) T r = T r (p) can be written for the 
case of Houston's experiments. In Houston’s experiment, P^g = P^. and T c . = T r . 
Heimel's empirical expression (eq. ( 8 ) in ref. 23; best fit to experimental data) that 
makes use of recent vapor pressure data can be approximated for P Cg in the range of 
Houston's data (area of thermionic interest) in the cesium arrival rate equation to give 


p = C exp 


(-9027/T ) 


so that 


in (-) 

\9027 / \C/ 


(6) 
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where C = 1. 193xl0 30 T”*"® 2 , T_ is the mean T in table XV:2 of reference 24 for the 

r r r 2 1 
range of arrival rates of interest, the units of p are particles (cm )(sec ), and T 

and T r are in K. Values of C for use in converter design are given in table HI, to- 
gether with the maximum error in ju caused by the use of T r rather than T,, in C. 
Combining equations (1), (2), and (6), we have 

cp= -/L^\ln + b * (7) 

\9027 ) \Cj 


and 


J = 120T 2 



1.2856 


a exp*' 11 605 b / T ) 


( 8 ) 


Equations (7) and (8), used in conjunction with equation (5) will be useful whether the 
cesium arrives from a gas under Knudsen or continuous flow conditions or even if it 
arrives from the substrate of an impregnated cathode. 

Equations (5) and (7) can be used in a practical converter configuration. Rufeh and 
Lieb (ref. 25) have reported on a variable spacing converter using a (110) oriented vapor- 
deposited tungsten emitter. They found the bare work function of the emitter to be 
4. 78 eV. Knudsen's flow condition is assumed. With C = 2.067X10 2,7 (cm -2 )(sec _1 ) 

(see table III) (In 1/C = 62. 90) and cp = 4. 78 eV, <p(ju) can be determined from equa- 
tions (5) and (7). The result is shown to be in quite good agreement with the data in fig- 
ure 3. 

The substrate work function can be obtained from cesiated emission data by solving 
equations (5) and (7) for the bare work function in terms T, ju, and cp: 


<P 0 = 


0. 158 cp - — In — - 0. 714 
9027 C 


0. 176 cp + 0. 121 


(9) 


Since cesium vapor pressure tables are normally available for reducing diode data, it is 
convenient to define an effective reservoir temperature T*; it is the temperature of a 
pool of cesium which has an evaporation rate p(T* = T r only when j u. = jLt Q ). Then 
equation (9) becomes 
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( 10 ) 


CP 0 = 


0. 158 cp + JL - 0.714 
T* 


0.176 cp + 0.121 


Equations (9) and (10) are expected to be accurate over the same domain as equations (5) 
and (7), that is, for 2. 3 < cp £ 3. 3 eV and perhaps for 4. 1 < cpg < 6. 5 eV. 

Equations (9) and (10) can be used as a check on experimental data, like the bare 
work function reported by Jacobson and Campbell (ref. 26), for example. They list a 
substrate work function of 4. 76 to 4. 79 measured in a vacuum diode. It is assumed that 
M based on continuum conditions applies since d/X > 10. Thus, we can adjust for 
arrival rate based on Mq by m = |/T r /T Mq. Setting T f = 620 K at T = 1916 K gives 
an adjusted j u = 4.6x10^ (cm“^)(sec"l). For this adjusted ju, T* = 595 K. When the 
value for <£> = 2.8 is used and inserted along with m and T r into equation (10), it is 
found that cpQ = 4. 8. Similar agreement was found with other combinations of d, T, and 


This check on diode data is dependent upon the test environments that exist in the 
cesiated and the vacuum work function studies. Since Cs is a good getter for oxygen, 
there could be different amounts of oxygen contamination on the electrode surfaces in 
the cesiated and the vacuum work function tests if the electrode material or the test en- 
vironment was altered. 

The effect of oxygen on cesiated diode performance can be quite strong, as was 
pointed out as early as 1931 by Villars and Langmuir (ref. 27). This effect appears to be 
beneficial to thermionic converter performance. As a result, it would be well if the 
emission properties of the multiple additive thermionic converter could be formulated by 

a simple relation. 

The analysis presented in this report can be used for oxygenated cesiated diodes by 
altering the bare work function <£>q. Kitrilakis , Lieb, Rufeh, and Van Someren (ref. 28) 
found that oxygen addition changed the bare -tungsten-surface work function from 4. 62 to 
values between 5. 02 and 5. 32. Figure (V-l) from reference 29 is used to compare pre- 
dictions of equations (1) and (5) and of the Rasor, Warner correlation in figure 4. Most 
of the data for a variety of temperatures (T) are contained within the two lines obtained 
by analytical estimate from equations (1) and (5). As the figure shows and as Kitrilakis 
noted, the early Rasor, Warner formulations fail at lower effective work functions. 


THE ANODE 


The temperature requirements on the anode (<1200 K) are much less stringent than 
on the cathode. This, of course, gives us much more latitude. 
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The anode generally has many more surface impurities than the cathode; so 
Houston's results are not necessarily pertinent to real converter anodes. Still it would 
improve our understanding if we could predict the minimum work function of "clean" 
cesiated anodes. 

Although the minimum work function ^ is 1. 5±0. 1 eV for many substrates 
(refs. 14 and 19) it is desirable to have a means of predicting <P m j n generally, since 
we have some freedom in choosing anode materials. Swanson and Strayer (ref. 13) have 
empirically determined that 


Acp min = 1.09(1.78 - cp Q ) (11) 

or 

<?0 + A( Pmin = <?min = L 9 5 ~ 09 <?0 C 12 ) 

Equation (11) was obtained by correlating field -emission data for Cs on molybdenum 
(Mo) , W, W(100), W(110), rhenium (Re), tantalum (Ta), and nickel (Ni). Although 
thermionic measurements are done at much lower fields and higher temperatures, 
table IV shows that Houston's results are in good agreement with equation (11). Thus, 
it appears that we may be able to use equation (12) for thermionic work functions also. 

Equation (12) can also be compared with experimental data. Jacobson and 
Campbell (ref. 27) have plotted their experimental data for <p m j n , and it is shown in 
figure 5 along with the prediction of equation (12). Fair agreement is obtained, Jacobson 
and Campbell data can be correlated by 

^min = 2 - 820 -°- 274 <?0 (I 3 ) 

Several investigators (refs. 29 to 31) have discussed the difficulties involved in de- 
termining collector work functions in thermionic converters. Hansen (ref. 30) suggested 
that the error band for the measured value of collector work function could be 0. 2 eV. 
Rufeh and Lieb (ref. 31) showed similar discrepancies in their measured values of col- 
lector work functions when compared with work function values actually achieved in 
practical regions of thermionic converter operations. 


DISCUSSION 

The accuracy of our formulation can be compared with that of the Levine (or Steiner) , 
Gyftopoulos and the Rasor, Warner correlations using Houston's data. The Steiner, 
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Levine, Gyftopoulos (SLG) correlation require knowledge of not only the bare work func- 
tion but also the monolayer work function and cesium coverage. As Houston did not 
measure these last two quantities, the SLG correlations unfortunately cannot be compared 
unambiguously. Houston (ref. 17) compared the calculated (Levine, Gyftopoulos) S- curve 

properties listed in table 1 of reference 6 with his data. (Table I was made up assuming 

14 2 

a monolayer work function of 1. 81 eV and a monolayer coverage of 4. 8x10 cm for Cs 
on 21 metals.) Houston concluded that "little correlation exists between the experimental 
results and the theoretical prediction. " 

The relations derived here do not include the effects of strong electric fields. In 
taking this into account Langmuir and Taylor (ref. 1) found that such fields affect emis- 
sion from cesium coated tungsten more strongly than that from pure tungsten - and 

variably with coverage. 


EXAMPLE APPLICATION OF THE FORMULATION 


Finally, to illustrate the use of the preceding equations the zero-field electron 
emission characteristics are calculated for a converter with the following properties: 

(1) Bare work function of the emitter, ^ = 4.9 eV 

(2) Bare work function of the collector, cp^ = 5. 0 eV 

(3) Emitter temperature, T g = 1900 K 

(4) Collector temperature, T„ = 850 K 

^ on 2 1 

(5) Cesium particle arrival rate at the emitter , f i = 1.43x10 (cm" )(sec ) 

From equations (5), a = 1. 42 eV and b = -1. 87 eV. Then from equations (7) and (8) 

we have (using C = 1. 794xl0 2 ^ (cm“ 2 )(sec“ 1 )) cp = 3.02 eV and J = 4. 35 (A) (cm -2 ), 
where cp and J are, of course, the cathode effective work function and electron 
emission current density, respectively. 

If we have a clean collector surface, equation (12) gives its minimum work function 


as - cp, 


mm 


1.49 eV. 


CONCLUDING REMARKS 

An empirical converter surface physics formulation has been provided which is 
quite simple and yet agrees better with experiment than other correlations in general use. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 9, 1971, 

120-27. 
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TABLE I. - RESULTS OF LEAST-MEAN-SQUARE FIT FOR a AND b 


FROM HOUSTON'S DATA (REF. 19) 


Material 

Constant, 

a, 

eV 

Constant, 

b, 

eV 

rms deviation, 

A rras’ 

eV 

Maximum deviation, 

A max J 

eV 

Range in work function 
<p compared, 
eV 

Pt 

1.13 

-1.401 

0.029 

0.044 

2.10 to 3.41 

Ir 

1. 33 

-1. 87 

.014 

.024 

2.04 to 3.45 

Os 

1. 35 

-1.83 

.014 

.033 

2.02 to 3. 59 

Rh 

1.42 

-2.09 

.022 

.052 

2. 31 to 3.60 

Ru 

1.43 

-1.93 

.031 

.076 

2.06 to 3. 55 

W 

1.48 

-1.88 

.013 

.026 

2. 14 to 3.26 


TABLE H. - ACCURACY EVALUATION FOR EQUATIONS (5) 


Material 

rms deviation, 

A rms’ 

eV 

Maximum deviation, 

A max’ 

eV 

Range in work function 
cp compared, 
eV 

Pt 

0.099 

0.17 

2. 10 to 3.41 

Ir 

.035 

.058 

2.04 to 3.45 

Os 

.016 

.028 

2.02 to 3. 59 

Rh 

.089 

. 17 

2. 31 to 3.60 

Ru 

.054 

.076 

2.06 to 3.55 

W 

.021 

.034 

2. 14 to 3.26 


TABLE m. - ARRIVAL RATE PARAMETER 


Cesium particle arrival rate, 

Arrival rate parameter, 

Maximum disagreement 

M, 

c, 

with HeimePs 

(cm~ 2 )(sec’‘ 1 ) 

(cm~ 2 )(sec“*) 

percent 

3.4X10 17 - 1. 1X10 19 

2.437X10 27 

8.9 

l.lxlO 19 - 1.4X10 20 

2.067X10 27 

7.6 

1.4X10 20 - 8.8X10 20 

1. 794X10 27 

6.6 
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Figure 4. - Experimental and theoretical 
cesiated-surface work functions with 
electronegative additive. 
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Figure 5. -Test of equation (12). 


E-6364 


NASA-Langley, 1971 22 


15 



National Aeronautics and Space Administration 
Washington, D. C. 20546 


OFFICIAL business 

PENALTY FOR PRIVATE USE $300 


FIRST CLASS MAIL 



POSTAGE AND FEES PAID 
NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 


POSTMASTER: 


If Undeliverable ( Section 158 
Postal Manual ) Do Not Return 


"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knoivl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concernmg its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 



TECHNICAL REPORTS: Scientific and 
echnicil information considered important, 
complete, and a lasting contribution to existing 

knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include conference proceedings, 
monographs, data compilations, handbooks, 
sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs,, 
Technology Utilization Reports and 
Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 


